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Abstract 

We give another derivation of quark confinement in QCD from the viewpoint of the 
low-energy effective Abelian gauge theory of QCD obtained via Abelian projection. 
It is based on the recently discovered Berezinskii-Kosterlitz-Thouless phase transition 
in the four-dimensional Abelian gauge theory. Moreover, we show that there exists a 
critical gauge coupling constant in QCD, above which confinement set in and below 
which there is no confinement. In a SU(N) gauge theory with Nf flavored fermions, 
we argue that this leads to a critical value of fermion flavors Nj for the confinement 
as well as the chiral symmetry breaking, which separates the deconfining and chiral 
symmetric phase from the confining and chiral symmetry breaking phase. A finite- 
temperature deconfinement transition is also discussed briefly. 
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Proving quark confinement in quantum chromodynamics (QCD) is a long-standing 
problem in theoretical particle physics. In a series of papers JD, 0, f5| [|, the author 
claimed that a novel reformulation of Yang-Mills (YM) theory leads to quark con- 
finement in QCD. The reformulation is based on the identification of the YM theory 
as a perturbative deformation of a novel topological quantum field theory (TQFT). 
The TQFT is suggested from the gauge-fixing and the Faddeev-Popov ghost terms in 
the maximal Abelian gauge (MAG) ||. The MAG is a choice of partial gauge fixing 
which realizes in the context of quantized gauge field theory an idea of Abelian projec- 
tion proposed by 't Hooft ||. After the Abelian projection, the original non- Abelian 
gauge theory is rewritten into the Abelian gauge theory with magnetic monopoles 
and adjoint matter fields. The TQFT in question is a field theory which describes 
the topologically non-trivial (background) gauge field configurations leading to non- 
vanishing magnetic monopole current loop (due to topological conservation law) in 
four space-time dimensions. As demonstrated in M, such topological field configura- 
tions in the functional integral give the dominant contribution to the string tension 
cr (or area decay of the Wilson loop), while the remaining part is identified with 
the perturbative contribution around such topological configurations (leading to at 
most perimeter decay). This result supports the dual superconductor picture of QCD 
vacuum JJJ. 

The term 'Abelian projection' is somewhat misleading, since it sounds as if off- 
diagonal parts are neglected a priori. Indeed, such a standpoint was taken in the 
early stage of investigation on quark confinement || [|. It is a trivial statement 
that the Abelian gauge theory is obtained from an non-Abelian gauge theory by 
neglecting the off-diagonal gauge fields. Note that the Abelian projection here takes 
fully into account the off-diagonal parts by integration (in the sense of Wilsonian 
renormalization group), which reduces to the renormalization of the effective Abelian 
gauge theory as shown in |1J. 

We adopt as a criterion of quark confinement the area law of the large (non self- 
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intersecting) Wilson loop (large compared with the typical size of QCD). What we 
have shown in the previous papers and what we are going to show in this paper 
are concerned with the physics in the low-energy or long-distance region of QCD. 
Therefore, we don't say anything novel about the high-energy or short-distance region 
of QCD where the asymptotic freedom is essential. This way of thought constitutes 
the core of Abelian dominance in low-energy physics of QCD ||. Actually, it turns 
out that Abelian dominance is lost in the high-energy or short-distance region of 
QCD. 

The above reformulation does not necessarily guarantee the non-existence of ad- 
ditional non-perturbative effects (coming from the remaining gluon self-interactions) 
which might invalidate the proof claimed above, although the validity of the identifi- 
cation is checked by numerical simulations. In this paper, in order to get rid of this 
doubt, we give another derivation of quark confinement by combining previous works 
and % § |. Q 

First, we consider a non- Abelian gauge theory described by the Yang-Mills (YM) 
theory with a compact (semi-simple) non-Abelian gauge group G. Then we decom- 
pose the gauge field into the diagonal and off-diagonal pieces (Cartan 
decomposition). The diagonal gauge field belongs to the maximal torus group H 
of G, while the off-diagonal gauge field to the coset G/H. The Abelian projection 
is realized as a partial gauge fixing which fixes the gauge degrees of freedom corre- 
sponding to G/H. Under the residual gauge transformation of H, the diagonal gluon 
fields transform as Abelian gauge fields, whereas the off-diagonal gluons transform as 
adjoint matter fields. Second, we obtain the Abelian gauge theory with a gauge group 
H by integrating out all the off-diagonal gauge fields belonging to G/H [|10|. This 

is possible by introducing the auxiliary antisymmetric tensor field B^ v . We called 

1 The main idea is still the dual superconductor picture of QCD vacuum where the magnetic 
monopole plays the crucial role. However, this does not imply that the instantons and other topo- 
logical non-trivial configurations such as vortex have nothing to do with quark confinement. Exact 
separation of monopole from other configuration is meaningless, since they are overlapped with each 
other. 
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the Abelian gauge theory obtained in this way the Abelian-projected effective gauge 
theory (APEGT) 0. 

A quite important observation is that the APEGT obtained in this way is consid- 
ered as the low-energy effective gauge theory (LEEGT) of the original non-Abelian 
gauge theory, e.g. QCD. This is because the off-diagonal gluons become massive after 
the MAG, as shown by direct numerical simulations |Tl| and an analytical consider- 
ation [Q]. Therefore, if the integration over the coset components performed above is 
identified with the integration over the high-energy mode in the sense of Wilsonian 
renormalization group, the resulting Abelian gauge theory is regarded as the LEEGT 
which is meaningful in the low energy- momentum region (p < O(m^)) or in the length 
scale £ larger than the inverse of the off-diagonal gluon mass {£ > 0{m~^) ). This 
is another important point which is sometimes misunderstood by the literatures. 

The APEGT of the YM theory with a gauge group G obtained according to the 
above procedures has the action |TJ (in 3+1 space-time dimensions), Q 



S APEGT — J 



d X 



l -^ L f»»F v - l -^g 2 B, v B^ + -z k B, v e, upa r + 



where f pu is the Abelian field strength f^ u := d^a u — d v a^ B pv is the rank 2 anti- 
symmetric tensor field and the dots denote ghost self-interaction terms and higher 
derivative terms 0. The renormalization factors are given by 

22 a 2 

z f = C 2 (G)Jl_lii//, (2) 

z B = +2C 2 (G)-f^ln/i, (3) 

z k = +4C 2 (G)-fVn/i, (4) 

where [i is the relevant scale of the theory and CiiG) is the quadratic Casimir op- 
erator, C 2 (G)5 AB = fACDjBCD_ The APEGT (D is an Abelian gauge theory with 
the gauge group H e x H m where H e = [/(l)^ -1 denotes the maximal torus subgroup 

of G (the subscript e is an abbreviation of electric) and H m is the same group as 

2 Although, in the paper [Q, the actual calculation has been presented only for G = SU(2), the 
extension to G = SU (N) is straightforward, at least in the one-loop level, as demonstrated in p2| . 



H e but identified with the gauge group for the dual gauge field (the subscript m 
is an abbreviation of magnetic). Here the dual gauge field 6 M is obtained from the 



Hodge decomposition of B^ v in 3+1 dimensions, B^ v — b^ u + x^vi 



Xfiu = \t l iva(}{p OL X^ — d a x 13 )- Consequently, the dual gauge field gets coupled with the 
monopole current as 



d A x B^r = / d*x B^r 



(fx b^ + 



(5) 



where the non- vanishing magnetic monopole current is deduced as violation of Bianchi 
identity for the Abelian part, := d u fuv, where /„„ := he^ upa f pa . 



YM Theory (G) 

Abelian projection 
APEGT (H e x H m ) 

J[db» 

Low-Energy 
Effective Theories 



Sym[A] 



SAPEGT[a,b,k] 




Integrating out off- 
diagonal parts G/H 



S m [k] 



Szsm [a] 



Monopole action 




Dual description 



SDGL[b] 



ZSM theory (H e ) Dual GL theory (H m ) 

The APEGT can give dual descriptions of the low-energy physics in the following 
sense. After eliminating a M field by integration, the APEGT yields the dual Ginzburg- 
Landau theory (DGL) with an action Sdgl^] which is written in terms of the dual 
gauge field 6 M alone. On the other hand, the APEGT is rewritten into another 
Abelian gauge theory with an action Szs[a] written in terms of the a M field alone after 
integrating out the b^ field. In this sense APEGT is an interpolating theory between 
two theories just described (We can also write the monopole current action S m [k]). 

The duality of the APEGT is explicit in the action ([!]), since at the level of field 

3 The theory was called the Zwanziger-Suzuki-Maedan (ZSM) theory, since it was obtained by 
Suzuki for G = SU(2) and Maedan-Suzuki for G = SU(3) |j] using the Zwanziger formalism |TJ| . 
However, the derivation of ZS theory was based on the ad hoc treatment of the off-diagonal parts 
under the assumption of Abelian dominance Q and assumption of the monopole condensation which 
leads to the mass for the dual gauge field as a result of spontaneous breakdown of the dual gauge 
symmetry H m (dual Meissner effect or dual Higgs mechanism) . These points are improved in jj| . 
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equation, B^ v = f^", or = B^, up to a multiplicative constant (renormalization 
factor) . 

It is possible to discuss the confinement through the DGL theory which is derived 
from APEGT as a LEEGT of QCD. However, in order to conclude quark confinement 
in this scenario (dual superconductivity f7|), it is necessary to show that monopole 
condensation really takes place. The analytical studies on quark confinement so far 
have been based on the assumption of monopole condensation in QCD. Although it 
is actually shown to occur at least in the strong coupling region based on the energy 
(action)-entropy argument in the lattice regularization [Q], we will adopt another path 
in this paper. In ([[]), g denotes the bare coupling constant and hence the prefactors 
Zf,ZB imply the renormalization effect. The calculated prefactor zb shows that the 
DGL theory Sdgl^ does not exhibit asymptotic freedom, as expected. 

On the other hand, eliminating the field B^ v from the interpolating theory ([!]) by 
a Gaussion integration, we obtain up to one-loop level Q 

which is nothing but the Abelian gauge theory with a gauge group H = [/(l)^ -1 . 
However, a remarkable difference from the true Abelian gauge theory is that the cou- 
pling constant g(fi) runs according to the same renormalization group beta function 
as the original YM theory [JJ (see eq. (fj)), 

/%) := = - Atf (A*) 3 , fro = UC2iG) > 0- (7) 

In this sense, this LEEGT reflects correctly the asymptotic freedom of the original 
YM theory. Indeed, the beta function (0) leads to the running coupling, 



4 Note that this action is valid just before the monopole condensation occurs. After the monopole 
condensation, the action should be modified. We consider approaching the confinement phase from 
the deconfinement phase. 
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As the length scale is larger and larger, the running gauge coupling becomes larger 
and eventually goes to infinity in the one-loop level, whereas, at two-loop level, it 

2 

converges to the infrared fixed point a = a* jjTJ| for a :— 

The purpose of this paper is to derive quark confinement in QCD based on 
the effective Abelian gauge theory (||) derived from the APEGT without assuming 
monopole condensation. We avoid to use the Zwanziger formalism which introduces 
an arbitrary unit four vector n^. Hence we do not need to worry about the ambiguity 
how to choose the vector n^. The same effect as monopole condensation is supplied 
with the vortex condensation based on the reformulation of the Abelian gauge theory 
(D as a deformation of another TQFT J3|. This approach can greatly simplifies the 
proof (e.g., we can avoid the tedious instanton calculations) and avoid the ambiguity 
just mentioned. In this strategy, we are dealing with the Abelian gauge theory alone. 
Hence the doubt on the identification of the perturbative sector does not occur, since 
the Maxwell theory has no self-interaction among the gauge fields. This is another 
advantage of this approach. Moreover, this view gives an insight into other interesting 
problems, e.g. quark deconfinement transition and chiral symmetry breaking (CSB) 



discussed below. Technical details will be given elsewhere [15 



Now we show that the quark confinement is realized in the non-Abelian gauge 
theory in the sense that the expectation value of the Wilson loop obeys the area law 
for sufficiently large (non self-intersecting) Wilson loop. In other words, the dominant 
contribution of the static quark potential between a quark and an anti-quark is given 
by the linear potential in the large separation (In addition, there is the Coulomb 
potential and a constant part which expresses the self-energy contribution). (] The 
information on the high-energy physics is incorporated into the APEGT through the 
running coupling constant (|§|). This is quite important to prove quark confinement 

based on the effective Abelian gauge theory. 

5 It should be remarked that we do not conclude anything new about the short distance behavior 
of the Wilson loop or static potential, because the following argument is based on the APEGT which 
is meaningful only in the large distance or low-energy region as explained above. 
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In fact, if we regard the above theory (H) with the LEEGT of the YM theory, it 
turns out that the quark confinement in QCD follows immediately from the result of 
the previous paper ||, i.e. the existence of a confining phase in QED. In || it was 
shown that the four-dimensional Abelian gauge theory has a confining phase in the 
strong coupling region g > g c where the fractional charge is confined by the linear 
static potential. The steps are as follows. The theory @ has the residual Abelian 
gauge invariance H. In order to quantize this theory, we must fix also the residual 
Abelian gauge degree of freedom. We adopt a manifestly covariant gauge fixing of 
Lorentz type. For a choice of gauge- fixing parameter (3 — — 2 0, the dimensional 
reduction occurs. Then the critical coupling g c is determined from the Berezinskii- 
Kosterlitz-Thouless (BKT) transition temperature Tbkt — 8 71 " of the two-dimensional 
0(2) nonlinear sigma model (NLSM): g c = tt. The weak coupling region g < g c is the 
conformal (Coulomb) phase. The above steps are shown schematically as follows. 

D-dim. QED 



Covariant gauge fixing 



D-dim. Perturbative QED 




deform 


D-dim. TQFT 


Dimensional reduction 


D-dim. Perturbative QED 


® 

deform 


(D-2)-dim. 0(2) NLSM 



Furthermore, the evaluation of the Wilson loop is performed using the equivalence 
of 0(2) NLSM with (neutral) Coulomb gas, sine-Gordon model and the massive 
Thirring model in two dimensions. This fact will shed light on the confining string 



picture [|16| in the low-energy region of QCD. 

Note that the Abelian gauge theory (|6|) is different from the usual perturbative 
QED in the following points. First, the Abelian group H is compact, since it was 
embedded in the compact (non- Abelian) group G, while the U(l) gauge group in the 
perturbative QED is considered to be non-compact. Due to compactness or peri- 



odicity, the Abelian gauge field can have the topologically non-trivial configuration 
(vortices). It is the periodicity why the fractional charge is confined while the integral 
charge is not confined. Second, the true Abelian gauge theory can not be asymptotic 
free. That is to say, in the infrared limit p j or long distance limit £ f <x>, the 
effective charge a(p) goes to zero. Therefore, QED inevitably goes into the Coulomb 
phase. On the contrary, the running gauge coupling (|§D of (§) increases monotoni- 
cally in the length scale £ and eventually becomes larger than the critical coupling 
a c with increasing £. Thus, in the low-energy region p < p c (or i > i c ) compatible 
with p < 0(mA) (or £ > (^(m^ 1 )), the theory (j6|) lies in the confining phase. The 
two critical scales m^, p c are expected to be the same order of magnitude (~ lGeV). 
Thus the critical scale p c where the confinement sets in is given by 

a{p c ) =a c = 7r/4. (9) 

The above consideration should be compared with the previous work || where 
quark confinement has been derived by directly studying the non-Abelian gauge the- 
ory without considering the LEEGT. The above consideration supplements the pre- 
vious treatment [0. By the existence of the scale p c , the separation of the gauge 
variable into and performed in has the following meaning. The back- 
ground piece fl^ represents the low-energy (or long-distance) mode p < p c while the 
perturbative piece V M the high-energy (or short distance) mode p > p c in QCD. 

From the result of || on the Abelian Wilson loop, therefore, the above considera- 
tion gives another proof of quark confinement in the sense of area law of the diagonal 
Wilson loop claimed in 0. Finally, the area law of the non-Abelian Wilson loop is 
shown as follows. The non-Abelian Wilson loop is rewritten as a path integral of the 
Abelian Wilson loop via a non-Abelian Stokes theorem as exemplified in the paper 
[0. Therefore, the area law of the non-Abelian Wilson loop follows from that of the 
Abelian-projected Wilson loop which is written in terms of the Abelian (or diagonal) 
part alone, see |T5| . 

In what follows, we consider a few applications of the above result. First, we 
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consider a SU(N) gauge theory with Nf flavored fermions. If the massless fermions 
are present in the above consideration, the Wilson loop loses its meaning as a criterion 
of quark confinement. This is because the virtual pairs of a fermion and an antifermion 
are produced from the vacuum and they break the string of inter-quark flux leading 
to the perimeter decay of the Wilson loop. Then the static potential is saturated to 
be a constant in the asymptotic region. However, if the fermion mass m/ is non-zero 
(even if small), the above criterion is still meaningful. The reason is as follows. Once 
the dynamical fermion mass m/ is produced, the fermions decouple below this scale, 
leaving the pure gauge theory behind. Therefore, the linear potential V(R) = aR 
survives the inclusion of light quarks for the length scale m^ 1 < £ < mfcr^ 1 (mj 1 a < 
p < niA). [] If the theory with fermions is asymptotic free and hence the running 
coupling constant is monotonically increasing toward the infrared region, a criterion 
of quark confinement is given by a(fi) > a c = | as long as the fermion is massive 

On the other hand, it is known that the chiral symmetry is spontaneously broken 
above the critical coupling a > a' c . For example, the SU(N) gauge theory has the 
critical coupling for the chiral symmetry breaking (CSB), a' c = 3C ^ R ^ = spf^ry 71 "? 
according to the naive Schwinger-Dyson (SD) equation [Q. This implies that the 
chiral symmetry is broken for not too large Nf below a critical value Nf, Nf < Nf. 
This leads to a critical value of fermion flavors NJ for the confinement as well as 
the CSB, which separates the deconfinement and chiral symmetric phase as follows. 
For N = 2,3, the critical coupling a! c of CSB is greater than or equal to the critical 
coupling a c of quark confinement, a' c > a c (In fact, for N = 2, a' c = An/ 9 and for 
N = 3 a' c = 7r/4 which happens to coincide with the critical coupling a c ). [] In the 
region a c < a* < a' c (a*: the infrared fixed point), although the fermion remains 

massless in the naive treatment without topological contribution, the vortex plasma 

6 Here we have implicitly assumed that the string tension does not change substantially despite 
the inclusion of dynamical fermions in the relevant region. 

7 In the region a' c < a < a c (which is expected to occur for large N) , the fermion is massive and 
the chiral symmetry is broken, but quark is not confined. 
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(monopole condensation) can induce the CSB which causes quark confinement. It is 
also known that monopole condensation enhances the CSB |]T7|. If this observation 



is correct, the CSB and confinement can occur simultaneously. 

We can estimate the critical number of flavors Nf for the deconfinement /confinement 
transition as follows. For SU(N) QCD with Nf fermions in the fundamental repre- 
sentation, the renormalization group beta function is given by 

dot 

(3(a) := fi— = -ba 2 (fi) - ca 3 ([i) ■ ■ ■ , (10) 

with the first two coefficients, 

1 1 / N 2 - 1 \ 

b = — (11N-2NA, c= 34JV 2 - l(WiV f - 3 N f . (11) 

6vr /; ' 24vr 2 \ N J K J 

Note that, above the value N f = (U/2)N {b < 0,c < 0) (N f = 11 for iV = 2 and 
Nf = 16| for N = 3), the theory is not asymptotic free, i.e., the beta function is 
positive for any a, (3(a) > 0. Below this value (6 > 0, c < 0), the infrared fixed 
point a* exists where the fixed point is given by a* = — b/c > 0. From this relation, 
we have a* = (^V^ty-^jv 3 ' so ^ na ^ ^ ne a b° ve relation is meaningful only when 
Nf > 34N 3 / (13N 2 — 3), since a* > only if this condition is satisfied (otherwise, 
a* < 0). Note that a* increases monotonically as the flavor Nf is decreasing and 
a* T +oo as N f | 34iV 3 /(13iV 2 - 3). For large N, N f approaches (34/13)JV. For 
Nf < 34iV 3 /(13iV 2 — 3), (3(a) < and the theory is asymptotic free (neglecting 
higher orders 0(a 5 )). This argument gives the lower bounds, Nj > 5.55 for N = 2 
and Nj > 8.05 for N = 3. The relationship between the critical flavor and the critical 
coupling for the deconfinement /confinement transition is obtained by putting a* = a c 
as 

= 2(22^ + 170^)^ (m 
f 8nN + a c (13N 2 -3)' 1 ' 

The critical coupling a c = tt/4 leads to Nf = 8.64 for iV = 2 and Nf = 11.9 for 
N = 3, under the assumption that the CSB has already occured. Therefore, the 
critical number of flavors is estimated: 8 > Nf > 6 for iV = 2 and 12 > Nf > 8 for 



N = 3. These results are consistent with the results of lattice simulations JIB . 
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Here we have used the result of SD gap equation for the CSB [jl4[. However, in 
order to really understand the relationship between confinement and chiral symmetry 
breaking, it is desirable to derive both of them on equal footing from QCD. This will 
be given in a forthcoming paper [15]. 

Next, we consider a finite-temperature deconfinement transition briefly. The ef- 
fective coupling at finite temperature behaves as o.eff{T) ~ l/(^o l°g(T/A)). If the 
temperature T is sufficiently high so that o.eff{T) < a c , quark confinement can 
not occur based on the above arguments. Due to monotonicity of qceff(T^ , there- 
fore, the confinement-deconfinement transition point is determined by the equation, 
aEFF{T) =a c = tt/4. 

The above strategy for proving quark confinement can be applied to the non- 
Abelian gauge theory defined in the D = d + 1 dimensional space time for arbitrary 
D > 2. For example, if we apply the above method to the 2+1 dimensional non- 
Abelian gauge theory, the 2+1 dimensional TQFT is , after the dimensional reduction, 
reduced to the 0+1 dimensional NLSM which is not a model of the quantum field 
theory, but that of quantum mechanics. It is nothing but the plane rotator model in 
quantum mechanics. In quantum mechanics, the phase transition does not occur and 
the theory lives in one phase. The phase should correspond to the confining phase 
of the 2+1 dimensional non-Abelian gauge theory. However, the Wilson loop can 
not be defined in 0+1 dimensional spacetime. We will need other criterion of quark 
confinement. 

Obviously, it is important to extend our argument to the two-loop level by deriving 



the APEGT of QCD at two- loop level [T5|, |I2" |. Furthermore, it is desirable to extend 
the arguments given for G = SU(2) in this paper to G = SU(N),N > 2. In a 
subsequent paper \TE\ we will deal with the generalized BKT transition corresponding 
to G = SU(N). Moreover, the area law of the non-Abelian Wilson loop will be derived 
based on fermionization of generalized sine-Gordon model (T5|. 

In H, quark confinement has been explained as a result of instantons in two- 
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dimensional NLSM. However, this is not inconsistent with the vortex picture for 
the effective Abelian gauge theory given in this paper, since it turns out that an 
instanton is composed of a vortex and an anti-vortex. The Abelian and non- Abelian 
quark confinement will be understood in a unified treatment [ |l5|] . 

This work is supported in part by the Grant-in-Aid for Scientific Research from 
the Ministry of Education, Science and Culture (10640249). 
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